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Abstract: An analytic method for predicting the performance of track-to-track association 
(TTTA) with biased data in multi-sensor multi-target tracking scenarios is proposed in 
this paper. The proposed method extends the existing results of the bias-free situation by 
accounting for the impact of sensor biases. Since little insight of the intrinsic relationship 
between scenario parameters and the performance of TTTA can be obtained by numerical 
simulations, the proposed analytic approach is a potential substitute for the costly Monte 
Carlo simulation method. Analytic expressions are developed for the global nearest neighbor 
(GNN) association algorithm in terms of correct association probability. The translational 
biases of sensors are incorporated in the expressions, which provide good insight into how 
the TTTA performance is affected by sensor biases, as well as other scenario parameters, 
including the target spatial density, the extraneous track density and the average association 
uncertainty error. To show the validity of the analytic predictions, we compare them with the 
simulation results, and the analytic predictions agree reasonably well with the simulations in 
a large range of normally anticipated scenario parameters. 

Keywords: track-to-track association (TTTA); sensor biases; analytic performance 
prediction; global nearest neighbor (GNN) 
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1. Introduction 

In distributed multi-sensor surveillance systems, the process of associating sets of local estimates 
from multiple sensors is a fundamental problem, named track-to-track association (TTTA) [1,2]. It 
is valuable for predicting the performance of TTTA, not only in multi-sensor tracking system design, 
but also in subsequent situation assessment and on-line resource allocation, such as sensor tasking. 
Generally, performance evaluation mainly depends on Monte Carlo simulations, which are expensive 
and time-consuming. Furthermore, little insight can be obtained for the intrinsic relationship between the 
relevant scenario parameters and the performance based on numerical simulations. Therefore, analytic 
prediction of data association performance, which does not resort to detailed and costly simulations, is 
of significant practical and theoretical importance. 

With respect to the analytic performance prediction of data association, Sea and Singer did pioneering 
work on the performance analysis of the nearest neighbor (NN) algorithm [3,4]. Saha provided a 
methodology for predicting the performance of a logic based ESM/radar TTTA algorithm [5]. Li et al. 
introduced the hybrid conditional averaging technique to predict the performance of the NN tracker [6] . 
Mei et al. gave the formula to calculate the theoretical probability of false association for TTTA 
using multiscan data [7]. Mori et al. derived an analytic expression, named the exponential law, to 
predict the performance of the global nearest neighbor (GNN) algorithm in terms of the probability 
of correct association [8,9]. On the basis of Mori's work [8], Ruan et al. developed an analytic 
performance prediction method for the feature-aided GNN algorithm [10]. Some oversimplified 
assumptions were employed in the derivations, including no false alarm, perfect detection, identical 
diagonal elements and all other identical elements in the feature confusion matrix, etc. Areta et al. 
presented procedures to calculate the probability that the measurement or the track originating from an 
extraneous target misassociated with a target of interest for the cases of NN and GNN association [11,12]. 
Mori et al. quantified the utility of additional feature/attribute information for TTTA in [13]. Although 
these works exhibit various degrees of success, they all do not take the impact of sensor biases on TTTA 
into consideration. 

In performing multi-target tracking with a sensor network, the systematic errors, i.e., sensor biases, 
can hinder the fusion of data from multiple sensors. Specifically, TTTA is always problematic with 
biased data, while sensor bias registration usually requires correctly associated data. That is, association 
and registration actually affect each other [14,15]. To fix this problem, many approaches were proposed 
to formulate association and registration as a joint optimization problem. However, these joint problems 
are very difficult to solve, due to the combinatorial nature of the problem and the non-convexity of the 
objective function [14-19]. Suboptimal solutions are obtained by performing association and registration 
in an alternately iterative manner. To a certain degree, the prediction of TTTA performance with biased 
data is also helpful for designing these joint association and registration algorithms. 

A number of TTTA algorithms are built upon the GNN criterion [1,2,16,20-23]. Most of them can be 
decomposed into two steps. The first step is to construct a matrix of association costs (or negative log- 
likelihood) between two sets of tracks. The second step usually consists of running a linear assignment 
algorithm on this cost matrix to determine the maximum likelihood association by globally minimizing 
the cost. The construction of the cost matrix was well discussed in [22-24]. For the second step, 
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many optimal linear assignment algorithms, such as Munkres [25], Jonker-Volgenant [26] and Jonker- 
Volgenant-Canstenon [27], can be adopted to identify the optimal associations. In this paper, we pay 
particular attention to predicting the performance of the GNN association algorithm with biased data 
analytically. It should be noted that Castella [28] did consider the influence of biases on the performance 
of the multiple site track correlation technique proposed by Singer and Kanyuck [29], but the TTTA 
method of Singer and Kanyuck applies a gating technique and makes the association decision based on 
the NN rule, rather than the GNN criterion. The goal of our study is to extend the results in [8,9], which 
established an effective approach for predicting the performance of the optimal assignment algorithms, 
in an attempt to incorporate the impact of sensor biases. 

The rest of the paper is structured as follows. Section 2 describes the single-scan TTTA problem. 
Section 3 presents the existing analytic prediction results of the bias-free situation. Section 4 provides 
the analytic prediction results for the GNN association algorithm with biased data. These results are 
verified via simulations in Section 5. Section 6 concludes the paper. 

2. Problem Statement 

In this paper, we consider a two-sensor single-scan TTTA problem. Although TTTA using multi-scan 
data has been theoretically proven to be better than that using single-scan data in [7], the former has 
higher computational complexity and may, in practice, exhibit lower power. That is because using multi- 
scan data increases the degree of freedom of the test statistics, which has a negative impact on the power 
of the test [30,31]. Therefore, we focus on the single-scan TTTA problem, which is a data association 
problem with two random point sets in an m-dimensional Euclidean space, lZ m . 

At a given time (without the time argument, for simplicity), the track set available at sensor A is 
denoted by (y is P^)^, where y< G K rn and P' A G n rnxm are the state estimate and error covariance 
matrix of the i-th target at sensor A, respectively; Na represents the number of tracks available at sensor 
A. Similarly to sensor A, the track set made by sensor B is (zj, P 3 B ) J^, where zj G TZ m and P B G 
TZ mxm are the state estimate and error covariance matrix of the j-th target at sensor B, respectively; 
Nb represents the number of tracks available at sensor B. Let y 0 and z 0 represent the dummy tracks of 
sensor A and B, respectively. When a track in a list is not associated with any track from the other list, 
it will be assigned to the dummy track. Following [8,32], we can define the association between the 
two track sets as a one-to-one mapping, a. Its domain is Dom(a) CI = {1, . . . , Na}, and its range is 
Rng(a) C J — {1, . . . , N B }. For i G Dom(a), a(i) is the index of the track from sensor B associated 
with track i from sensor A. Tracks in X \ Dom(a) are associated with z 0 , while tracks in J \ Range(a) 
are assigned to y 0 . 

Let a* be the ground-truth association mapping. The track model is given as follows. We assume that 
sensors have translational biases on their local estimates of targets. The biases are modeled as additive 
constants. For a common target, i, we have: 

yi = x» + b A + Ui, (1) 
z 0 . w = Xi + b B + Vi, (2) 
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where Xj is the real state of target i, b^, b B represent the biases of sensor A and B, respectively, u { , Vj 
are the random state estimation errors of target i at sensor A and B, respectively, and u« ~ A/"(0, P A ), 
v f ~JV(0,P&). 

The association can also be defined as an (N A + 1) x (Nb + 1) matrix h between (y^jJg an( ^ ( z i)j=0' 
with (i, j) element for z = 1, . . . , Na + 1, j = 1, ■ ■ ■ , Afe + 1, such that 

z/ y i( Zj are from the same target ^ 
0, otherwise 

The relationship between a and h is: if a(i) = j, then = 1; if % e X \ Dom(a), then ft, i Q = 1; if 
j E J \ Range(a), then ft 0 j = 1; otherwise ft^j = 0. 

The GNN algorithm determines the best association decision from all possibilities by minimizing the 
overall cost [1,2]: 

N A N b 




EE^- (4) 



mm 

h 

i=Q j=0 



subject to the constraints: 



N A 

E^i = 1- J/O (5) 

i=Q 

E^i = 1, (6) 

j=0 

hij E {0,1}. (7) 

Constraints (5) and (6) guarantee that each track from a list is assigned to one and at most one track from 
the other list or to the dummy track. In problem (4), Dy (i, j ^ 0) is the cost between y^ and Zj, defined 
by their normalized distance as: 

Dij = \\zj - yiHjLi, (8) 

where ||x|| M denotes a semi-norm of a vector, x, defined by a nonnegative symmetric matrix, M, as 
||x|| M = Vx T Mx; Sjj is the average association uncertainty covariance matrix and S» = + P^ — 
P AB ~ ^*ab '■> ^ab i s me cross-covariance matrix, representing the correlation in the track errors, which 
results from the common process noise [33]. Because the correlation between tracks is difficult to obtain 
in practice, one can ignore it by assuming P A 3 B = 0 [23,24]. The effect of the cross-covariances on TTTA 
is beyond the scope of this paper, which can be found in [33,34]. We further assume that P^ = P^ and 
P J B = P B for all i and j, such that S = P^ + P B is strict positive definite. 

In the GNN algorithm, the probability of correct association, Pc, of target i is defined as: 



c 



P{a(i) = a*(i)}, (9) 



where a is the GNN optimal association estimate. 
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3. Previous Results with Bias-Free Data 

Mori et al. [8,9] derived the analytic expression for the correct association probability, Pc, of the 
GNN algorithm with bias-free data (i.e., h A = 0 in Equation (1), h B = 0 in Equation (2)), which obeys 
an exponential law: 

P c = exp(-C m p T a m ), (10) 

under the assumption of no false alarm and perfect detection. Considering the effects of extraneous 
objects, the expression of Pc was changed into: 

Pc = exp[-{C m (3 T + D m (3 F )a m }. (11) 

The relevant parameters and assumptions used in Equations (10) and (11) (also effective in the 
subsequent analysis) are given as follows. 

• C m is a constant defined by: 

C m B m C m , (12) 

where: 

C m = 2 m - 1 vr- 1 /2r(^±A), (13) 
and B m is the volume of the unit m-dimensional ball: 

m 

7T 2 

Bm = r(f + 1)' (14) 

• (3t denotes the spatial density of targets defined as the expected number of targets in a unit volume 
of the measurement space, lZ m : 

B m r m 

Targets are assumed to be independently identically distributed with a common distribution, which 
is uniform on a m-dimensional ball, with a large enough radius, r. The target number, N T , is a 
Poisson random variable with a mean, \t, that is: 

p[NT) = e " T W: (16) 



D m is a constant given by: 



where: 



D m = B m D m} (17) 

Dm = 2 m/2 ^v- (18) 
1 {m/2) 

f3p denotes the spatial density of extraneous objects: 

Pf = — ■ (19) 

J - J m' 

The states of the extraneous objects share the same distribution with targets. The number of false 
objects, N F , is another Poisson random variable with mean X F , that is: 

P(iV F ) = e-^-f-. (20) 
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• a is the average innovation standard deviation. 

Actually, the results of Equations (10) and (11) are specified for measurement-to-track association in 
single-sensor multi-target tracking problems. The novelty of this paper is to extend these results to the 
TTTA problem with biased data in multi- sensor multi-target tracking applications. 

4. Analytic Performance Prediction of the GNN Algorithm with Biased Data 

In this section, we consider the impact of sensor biases and propose the analytic prediction results for 
the GNN algorithm. 

4.1. Association without Extraneous Tracks 

First, we start with the two-target case, target i and j (i ^ j). An optimal association, a, is correct, 
i.e., a(i) = a*(i) and a(j) = a*(j), if and only if the cost of the correct association is less than that of 
the incorrect association, that is: 

A,a*(i) + Dj >a *(j) < A,a*(j) + -Dj,a*(j), (21) 

or, equivalently, if and only if: 

AJy = ||z a *(j) — y»|||-i + ||z a *0') ~~ y i 1 1 s - 1 

i,a*(i) j,a*(j) 

-K*(j) - y*lls-i , - IK*« - y^lls-i , ( 22 ) 

= -2(z a . (i) - z a . U) )S-\yi - yj) < 0. (23) 

Considering track model (1) and (2), we have: 

z a *(i) - z a * (i ) = (xj + b B + Vj) - (xj +h B + Vj) = (xi -Xj) + (v« - v^), (24) 
yi-Yi = ( x i + b A + Ui) - (xj + b A + uj) = (xi - Xj) + (Uj - u_j). (25) 

Then, Equation (23) is equivalent to: 

(Aug - Av ii ) r S- 1 Ay y < HAy^Hl^, (26) 

where: 

Au^ = Ui- uj, (27) 

Avy = Vi-Vj-, (28) 

Ay,,- = yi-yj. (29) 

Thus, based on the detailed derivations in Appendix A, the probability of correct association with two 
targets is: 

7) ' ( 3 °) 
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where o is the average association standard deviation, and: 

a = [det(S)]3S. (31) 

It is very difficult to extend Equation (30) to the multi-target case, because the incorrect associations 
may have more complicated errors besides the two-object transposition, such as the multi-object 
transposition. To overcome this difficulty, let the event, Ec{a(i) = a*(i)} (target i is correctly 
associated), be approximated by the event that there is no two-object transposition of A > 0 involving 
target i. With Equation (16) and a further assumption that the event { A Jy > 0} is independent for each 
j, i.e., each "potential" transposition is independent, we have: 



oc 



<7 



P{Ki) = a*(i)\N T }P(N T ) (32) 

N T =0 
oo Nt 

II P{^<0}P(N T ) (33) 

N T =0j=l,jjH 
oo _ 

N T =0 



[l-C m (-) m f T ^P(N T ) (34) 



r 



where: 



and: 



, M N T =0l L Mfj J jV T ! 
l"^ m (f) m 

C , ^ n exp(-C' m /3 T ) (36) 
C' m exp(-C m f3 T a m ), (37) 



c: " i rUr : (38) 



(3 T = B m p T a m (39) 



is the normalized target density, which is defined as the expected number of targets in the one-sigma 
ellipsoidal volume derived from the average association uncertainty covariance matrix, S. The formula: 



-oo 

X 



8=0 

is used in the derivation of Equation (36). The idea of independent transpositions was first introduced 
in [35], which gave some justification of this assumption. Generally, all possible association mappings 
are considered to be independent with each other and totally random. In this sense, the independent 
transposition assumption is reasonable. 
Generally, r ^> a. Hence: 

C m w 1, (41) 
which was used in the derivation of Equation (10) in [8]. Consequently, Equation (37) becomes 

P c w exp(-C m p T a m ). (42) 
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To maintain consistency with [8-10], Equation (42) is used in the following discussions. Since the 
bias terms are canceled out in Equations (24) and (25), Equation (42) is the same as (10), which does not 
account for sensor biases. Then, we can see that without extraneous and missed tracks, the translational 
biases do not affect the performance of the GNN association algorithm. 

4.2. Extraneous Tracks Only Existing in One Track Set 

In Section 4.1, we assumed that the cardinalities of the two given random point sets are identical. In 
reality, false tracks may disturb the situation. In this section, we discuss the effects of false tracks on the 
TTTA performance with biased data, in which there are false tracks only in the track set of sensor B. The 
event, E c {a(i) = a*(i)}, can be written as 

E c {a(i) = a*(i)}^Ec T f)Ec F , (43) 

where Ec T is the event that there is no transposition between target i and any other target and Eq f is the 
event that target track yi is not associated with any false tracks. Since the targets and false alarms are 
independent from each other, we have: 

P c = P(E c {a{i) = a*(i)}) « P(E Ct )P(E Cf ). (44) 

We can use the result of Equation (42) for P(E Ct ). The focus now is to analyze P(E Cf ) in the situation 
where sensors have translational biases. 

Given an extraneous track, z e , and the target track, z a *^, from sensor B, let P ie be the probability that 
z e is not associated with yi , that is: 

P ie = P{||z a * (i) - yi |||_ a < ||z e - YilH-i}. (45) 

Without loss of generality, let y, ; = 0. Then, Equation (45) is equivalent with: 

P ie = PQz^Wl^ < Kill-!}. (46) 

According to the track model (1), we can get the real state of target i, for the form only, as: 

Xi = -b A - u;. (47) 

Consequently, together with Equations (2) and (47), we have: 

z a *(i) = -bA - Ui + be + Vj = Ab + Vj - u i; (48) 

where: 

Ab = b B - b A , (49) 

is the relative translational bias vector between the two sensors. Then, based on Equations (1), (2) 
and (48), we have z tt *(j) ~ A/"(Ab, S). Thus, ||z a *(j)|||_i is distributed according to the noncentral 
chi-squared distribution with the degree of freedom, m, and the noncentrality parameter, 5, i.e., 
H z a*(i)Hs-i ~ Xm(^) [36]. The noncentrality parameter, 5, is: 

5 = Ab^S^Ab. (50) 
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Assuming that the extraneous track, z e , is uniformly distributed on the m-dimensional ball with radius 
r, the distribution of ||z e |||_i can be expressed as: 

P{\W\\U<y} = {^ L ) m . (51) 

(7 

Let Y = ||z e |||_i; then, its probability density function is 

My) - - S&r""'. (52) 



dy 



Let X = ||z a .(j) |||_! and the probability density function of X be denoted by fx(x; m, 5). Considering 
that the extraneous track, z e , and the target track, z a *u\, are independent, the joint probability density 
function of random variables, X and Y, can be calculated as: 

fx,v{x,y) = fx(x)f Y (y) = f x (x;m,5)f Y (y). 

Then, we have: 

Pie = P{X<Y} (53) 

fx,y(x,y)dxdy (54) 




{{x,y):x<y} 

+00 f+OO 

fx(x)f Y (y)dydx (55) 



-OO J X 



m/2 



fx(x-m,5)^-y m l 2 - 1 dydx (56) 



m/2 r + ™ r^ )2 

- 1 f x (x;m,S) tf^dydx (57) 



JL\m 

a> J -co 

/ fx{x;m,5)^ dx (58) 

fx(x;m,5)(-) m dx- x m/2 f x (x;m, S)dx] (59) 



r \m 



1 



+oo 



I -TP- \ x m / 2 f x (x ]m ,6)dx (60) 



^ ~ 7^^l(s)(x m/2 ) (61) 



(62) 



where: 



D' m = £ x?n (8)(x m/2 ) = / x m ' 2 f x (x- m, S)dx. (63) 

J — oo 

If m = 2fc and A; is a positive integer, Equation (63) becomes 

D' m = £ xUS) (x k ), (64) 

which is the kth order moment of the noncentral chi square distribution with degree of freedom m and 
noncentral parameter 5. If m — 2k + 1 and A; is a nonnegative integer, Equation (63) is the fractional 
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moment of the noncentral chi square distribution, and the calculation method of Equation (63) can be 
found in [37]. 

When there are N F false tracks in the track set of sensor B, with the independence assumption, the 
probability that no extraneous track is associated with y, is 

P(Eo r \N F )=P»*. (65) 

Under the assumption of the false track number in Equation (20), we have: 

oo . n f 

P(E Cf ) = £ fjPZ' (66) 

N F =0 F ' 

= exp[-X F (l - P te )] (67) 
= exp[-\ F D m (°) m ] (68) 
= exp[-p F D' m } (69) 
= exp[-B m p F d m D' rn \ (70) 
= exp[-D' m P F a m ], (71) 

where 

0f = B m (3 F a m (72) 
is the normalized density of false tracks, and 

D' m = B m D' m . (73) 

By substituting Equations (42) and (71) into Equation (44), we have: 

P c = exp(-C m p T a m )exp[-D' m (3 F a m ] (74) 
= exp{-[Cj T + D'J F )} (75) 
= exp{-[C m (3 T + D' m P F ]a m }. (76) 

When there are no extraneous tracks, i.e., (3 F = 0, Equation (76) is the same with Equation (42). The 
effect of sensor biases is coupled with false tracks by the multiplying factor, D' m (or D m ). With m = 2, 
we have C m = tt, D' m = tc(2 + 5); then: 

P c = exp{ - tt[/3 t + p F (2 + 5)]a 2 }. (77) 

In addition, if sensors do not have biases (or have the same translational biases), i.e., Ab = 0, then 

5 = 0, D' rn = 2n, and 

P c = exp{ -tt[(3 t + 2(3 F }a 2 }, (78) 
which is consistent with Equation (11) under the two-dimensional case. 
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4.3. Extraneous Tracks Existing in Both Track Sets 

When false tracks exist in both track sets, with a simple generalization of Equation (43), we have: 

E c {a(i) = a*(t)}^E CT f]Ec FA f]Ec FB , (79) 

where Eq fa is the event that track y, is not associated with any false track from sensor B and Ec FB is 
the event that track Zj is not associated with any false track from sensor A. With the assumption that 
these events are independent, Equation (44) then is extended to: 

P c « P(E Ct )P(E Cfa )P(E Cfb ) (80) 
= exp{-[C m 0T + D' m (0 FA + p FB )]} (81) 
= exp{-[C m /3 T + D' m ((3 FA + p FB )}a m }, (82) 

where (5 FA and /3 FB are the densities of false tracks made by sensor A and B, respectively; (3 FA and (3 FB 
are the normalized densities of false tracks made by sensor A and B, respectively. 

Comparing Equation (82) with (76), it seems that the situation with false tracks in both track sets has 
the same difficulty with that in only one sensor. However, this is not true. As shown in Section 5.3, 
the relative quantities of false tracks between sensors have a large impact on the Pc\ so, the simple 
sum of (3 FA + (3 FB in Equation (82) is not sufficient to reflect the intrinsic impact of false tracks on the 
correct association probability. Other risks in analytic performance prediction using Equation (82) will 
be examined through simulations in Section 5.3. Here, we just give the main conclusion that in a large 
range of normally anticipated operating conditions, Equation (82) can be used as an upper bound for the 
performance of the GNN association algorithm in terms of the correct association probability. 

5. Simulation 

To show the validity of our theoretical analysis, we test the above analytic results by comparing them 
to the simulation results. We consider the single-scan TTTA problem and make use of the Munkres 
algorithm to solve the GNN assignment problem [25]. The scenario consists of two sensors tracking 
multiple targets. The targets are distributed uniformly in a two-dimensional ball space (m = 2) with 
target extent radius r. The number of targets is a Poisson variable with mean At = 20, and the 
target spatial density is varied by changing r. The random error covariance matrices of the target state 
estimates are presumed to satisfy = P B . The correlation among tracks from different sensors due 
to the common process noise is neglected in the following simulations, and the detailed analysis and 
simulations for the correlated tracks can be found in [33,34]. Sensors may report extraneous tracks, 
which are assumed to have the same spatial distribution with the target state. The track detection 
probability of each sensor is presumed to be one. Since the translational biases are in Cartesian 
coordinates, only the difference of them is observable. Therefore, the relative biases between sensors, 
instead of the absolute biases of each sensor, are considered in simulation. Let Ab x , Ab y represent the 
relative translational biases in the x,y coordinate, respectively. All of the simulation results are based 
upon 500 Monte Carlo runs. 
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5.1. No Extraneous Tracks 

When we study the validity of Equation (42) with biased data and no extraneous tracks, two simulation 
experiments are designed. The target spatial density, /3t, is set to be 0.01, 0.025, 0.05 reports /km 2 , with 
the target extent radius r being 25.23, 15.95, 11.28 km, respectively. In the first experiment, the average 
association uncertainty covariance matrix is fixed to be S = diag[l, 1] km 2 , while Ab x , Ab y vary from 
—3 km to 3 km in increments of 1 km. Tables 1, 2 and 3 show that the correct association probability, 
Pc, obtained by simulations remains approximately unchanged and matches the analytic prediction well, 
even if the sensor biases vary in a big range. That is because, when no extraneous tracks exist, the Pc 
of the GNN association algorithm is not associated with translational biases in Equation (42). From the 
comparison between Tables 1, 2 and 3, the theoretical Pc value given by Equation (42) is less accurate 
as the track density increases. An explanation of this phenomenon is that, although the assumption of 
a large enough r is necessary in the derivation of Equation (42), the target extent radius, r, is not large 
enough under a high target density. Another possible reason might originate from the more complicated 
association errors, i.e., multi-object transpositions, which we did not model. However, the statistics 
concerning such "high-order" transpositions are not clear at the moment. In the second experiment, 
sensor biases are fixed to be Ab x = Ab y = 1 km, while the average association variance, a 2 , varies 
from 0 km 2 to 2 km 2 in increments of 0.4 km 2 . Figure 1 plots the correct association probability, 
Pc, as a function of a 2 , and each point is compared with the prediction (solid line) by Equation (42). 
The analytic prediction results match the simulation results very well, with a few exceptions: when the 
targets are dense and the average association variance is large. In view of the approximations we have 
employed, including some quite "bold" ones, such as the independent transposition assumption, this 
consistency between the simulation results and the theoretical prediction is rather remarkable. 

Figure 1. Correct association probability as a function of the average association uncertainty 
variance — analytic prediction results against simulation results. Target spatial density ft T is 
0.01, 0.025, 0.05 reports /km 2 . Sensor biases are set as Ab x = Ab y = 1 km. 
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Table 1. The impact of sensor biases on correct association probability — the prediction by 
Equation (42) is 96.90%. (3 T = 0.01 reports /km 2 . 



Pc (%) 


— OK.II1 


9L-m 


— Ikiu 


Ab x 
0km 


1 K 1 1 1 
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Ab y 0 km 
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1 km 
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97.20 
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97.62 


97.67 


97.32 


3 km 


97.26 


97.08 


97.28 
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97.19 


97.17 


97.11 



Table 2. The impact of sensor biases on correct association probability — the prediction by 
Equation (42) is 92.44%. (3 T = 0.025 reports /km 2 . 
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Table 3. The impact of sensor biases on the correct association probability — the prediction 
by Equation (42) is 85.46%. (3 T = 0.05 reports /km 2 . 
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5.2. Extraneous Tracks Only in One Track Set 

When we test Equation (76) with extraneous tracks in the track set of sensor B, the normalized target 
density, $ T , is fixed to be 0.05 reports, with the average association uncertainty covariance matrix, 
S = diag[0.5, 0.5] km 2 , and the target extent radius, r = 14.14 km. The relative translational biases 
are set as Ab x = Ab y = 0, 1, 2 km. The ratio, $f/$t, is changed from zero to one in increments of 
0.2. In each Monte Carlo run, false tracks are added according to given densities. Figure 2 shows the 
correct association probability as a function of the ratio, $f/$t, which determines the relative density 
of extraneous tracks to target tracks. The analytic prediction results (solid lines) match the simulation 
results (points) very well, and they perform better matching when small sensor biases occur. 

Figure 2. Correct association probability as a function of ftp / fix — analytic prediction results 
against simulation results. The normalized target spatial density, (3t, is fixed to be 0.05 
reports. Sensor biases are Ab x = Ab y = 0,1,2 km. 
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5.3. Extraneous Tracks in Both Track Sets 



When we test Equation (82) with extraneous tracks in both track sets, the average association 
covariance matrix is set to be S = diag[0.5, 0.5] km 2 , and Ab x is assumed to be equal with Ab y varying 
from -3 km to 3 km in increments of 1 km. Let (3p be the total normalized density of extraneous tracks, 
and /3p = (3fa + $fb- Let a be the ratio of false tracks from sensor A, and a = $fa/$f- Given the total 
density of false tracks, the parameter, a, reflects the allocation (or the relative quantity) of the extraneous 
tracks between sensor A and B. Four cases of a (a = 0, 0.1, 0.3, 0.5) are considered in the simulations. 
To better understand the prediction performance of Equation (82), two groups of simulation experiments 
are performed. In the first group, the normalized target spatial density, p T , is fixed to be 0.01 reports 
(accordingly, r is set as 31.62 km). Figures 3 and 4 plot the correct association probability as a function 
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of sensor biases, with f3 F = 0.1/3 T and /3 F = f3 T , respectively. In the second group, f3 T is fixed to be 
0.05 reports (accordingly, r is 14.14 km). Figures 5 and 6 plot the correct association probability as a 
function of sensor biases, with (3 F = 0.1(3t and $ F = (3 T , respectively. 

Figure 3. Correct association probability as a function of sensor biases when both track 
sets contain extraneous tracks — analytic prediction results versus simulation results. The 
normalized target spatial density, $ T , is fixed to be 0.01 reports. The total normalized false 
track density is set as (3 F = 0.1/3r. 
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Figure 4. Correct association probability as a function of sensor biases when both track sets 
contain extraneous tracks — theoretical prediction versus simulation results. The normalized 
target spatial density, f3 T , is fixed to be 0.01 reports. The total normalized false track density 

is set as f3 F = /3 F - 
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Figure 5. Correct association probability as a function of sensor biases when both track 
sets contain extraneous tracks — analytic prediction results versus simulation results. The 
normalized target spatial density, $t, is fixed to be 0.05 reports. The total normalized false 
track density is set as (3 F = 0.1(3 T . 
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Figure 6. Correct association probability as a function of sensor biases when both track 
sets contain extraneous tracks — analytic prediction results versus simulation results. The 
normalized target spatial density, j3 T , is fixed to be 0.05 reports. The total normalized false 
track density is set as ftp = $t- 
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Since a series of assumptions are made in the derivation of Equation (82), the prediction accuracy 
degrades in the scenario where false tracks exist in both sensors. When the number of extraneous objects 
is far less than the number of targets 0 F = Q.1(3 T ), the prediction results are an upper bound for the 
simulation results, and the prediction deviation is below 7%, as shown in Figures 3 and 5. While there 
are a lot of false tracks ({3 F = (3t) in Figures 4 and 6, Figure 4 shows that, in the sparse-target scenario 
(Pt = 0.01 reports ), even with lots of false tracks, Equation (82) can also be used as an upper bound 
for the simulation results. Comparing Figure 6 with Figures 3-5, we can see that this prediction upper 
bound becomes invalid only in the scenario with a high target density (J3t = 0.05 reports), a high false 
track density (/3 F — Pt) an d large sensor biases (| Ab x \ > 1.3km). Meanwhile, the prediction departure 
becomes larger as a increases. The fundamental reason for the prediction departure is not clear at this 
point. Two possible explanations are given as follows. First, in this complicated scenario, especially 
with a large a (< 0.5), false tracks of different sensors may associate with each other. This leads to a 
"reduction" of the false track density. In this sense, fi F A +Pfb in Equation (82) overestimates the density 
of false tracks. In addition, this false track "reduction" effect can be amplified under large sensor biases, 
since the multiplying factor, D' m , in Equation (82) is a quadratic function of sensor biases. Second, when 
both track sets have false tracks, the multi-object transpositions we did not model are apt to occur. Yet, 
for all that, under large sensor biases and a large a (in Figure 6), the extreme pessimistic tendency of the 
prediction does not matter very much, because it appears only when Pq is quite small (below 50%). 

To sum up, in a large range of normally anticipated operating conditions, Equation (82) can be used as 
an upper bound for the simulation results. Further investigations are needed to provide a more accurate 
prediction of TTTA in the cases where false tracks exist in both sensors. 

6. Conclusion 

In this paper, we proposed an analytic performance prediction method for the GNN association 
algorithm with biased data in multi-sensor multi-target tracking applications. The probability of correct 
association is adopted as the performance criteria. The novelty of the method lies in that it accounts for 
the sensor biases in the TTTA problem, and analytic approaches are developed to reveal the intrinsic 
relationship between a set of key scenario parameters and the performance of the optimal TTTA 
algorithm. To verify the validity of the predictions, we compared them with the results of Monte-Carlo 
simulations. These shows that the analytic predictions agree reasonably well with the simulation results. 
In a word, the main results of the paper are as follows: 

(1) When both track sets are free from extraneous tracks, the correct association probability of the 
GNN association algorithm in Equation (42) has nothing to do with the translational biases of sensors. 
It is only related with target spatial density and the average association errors. 

(2) When one track set suffers from extraneous tracks, an exponential law was given in Equation (76), 
which uncovers the relationship between the correct association probability and the scenario parameters, 
including the average association errors, sensor biases, the target spatial density, as well as the false track 
density. The impact of sensor biases is embodied by an amplification coefficient of the false track density. 
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(3) When both track sets contain extraneous tracks, an analytic upper bound of the correct association 
probability for the GNN association algorithm was proposed in Equation (82), which is effective in a 
large range of scenario parameters. 

A series of assumptions were made to get concise and explicit mathematical results, including 
no missed detection, the translational bias model, identical association uncertainty covariance on all 
tracks, etc. The results derived based on these assumptions are not intended to be used directly in 
practice. The major contribution of the paper is to help designers understand the intrinsic relationship 
between a set of key scenario parameters and the performance of TTTA and to demonstrate the 
possibility that analytic performance prediction can be a potential substitute for the costly Monte Carlo 
simulation method. 

To make the proposed analytic method more practical, a lot of work is needed in the 
following directions: 

(1) Further investigations are needed to improve the prediction accuracy of our method in the scenario 
where false tracks exist in both sensors. The statistics concerning the more complicated association 
error,s such as multi-object transpositions, can be addressed in the future. 

(2) For simplicity, we assumed that the track detection probabilities of both sensors are one. However, 
missed detections and false tracks do happen in real environments. To complicate matters, when missed 
detections occur at one sensor, the corresponding tracks that are supposed to associate with these missed 
tracks in the track set of the other sensor are equivalent to extraneous tracks for the remaining tracks of 
the former sensor. When sensors suffer from both missed detections and extraneous tracks, accounting 
for the missed detections in the formula of the correct association probability needs to be investigated in 
the future. 

(3) We only examined the impact of translational biases on the performance of TTTA. Since the 
impact of translational biases on the tracks is uncoupled with the states of targets, the translational 
biases can be named state-independent biases. However, in some applications, sensors can have biases 
on their range and azimuth measurements. The impact of these biases on the position estimates of 
targets is dependent on the real states of targets. Therefore, we name the range and azimuth biases 
the state-dependent biases. Under the scenario with state-dependent biases, the bias terms cannot be 
canceled out in Equations (24) and (25). Thus, the state-dependent biases affect the performance of the 
GNN association algorithm, even if there are no false tracks and missed detections. How to deal with 
the more complicated state-dependent biases in analytic performance prediction is left for the future. 
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Appendix A 



In this appendix, we give the derivation of Equation (30). Most of the results in this appendix are 
taken from [8]. 

P c = P{a(i) = a*(i)} = P{AJ ij <0} (83) 
= £[P((A Uij - Av u -) T S- 1 Ay JJ < ||Ay y |||_ 1 )|Ay^] (84) 
= 8(erf(a i:j )) (85) 

where = || Ay^ ; || s -i/v2, £[■} is the mathematical expectation operator, erf(-) is the error 
function defined as erf(c) = ^- f° e~ t% l 2 dt and its complement is denoted by erfc(c), i.e., 
erfc(c) = 1 — erf(c). 

Consider the case where target i is located at the origin, while target j is at a random position 
according to the uniform distribution within the m-ball with radius r. Then, we have: 



£(erf( aij )\ yi = 0) = 1 - {B^Y 1 ! er fc( aij )d Ay tJ (86) 

= 1 - (B m r m )- 1 2 m / 2 [det(S)] 1 / 2 / erfc(\\ V \\)drj (87) 

JV{r) 

« 1 - B-Wi-r [ erfcdlvlDdrj, (88) 

r JTZ m 



where V{r) = {t] G TZ m \^/2r] T Sr] < r}; the approximation was made by replacing V{r) with TZ m , 
assuming that r is large enough. By using a kind of spherical integral shown in Appendix B, we have: 

/* /*oo 

/ erfc{\\ V \\)dr] = mB m / p^erfcpdp (89) 
Jn m Jo 

= B rn — / p m e- p2 ' 2 dp (90) 
2n In 



B^I^VC^). (91) 



Then, based on Equations (88) and (91), we have: 



S(erf (a^lyi = 0) = 1 - C m (-) m , (92) 

r 
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where C m was given in Equation (13). Furthermore, since r is large enough, we can replace the 
condition, y, = 0, by "when the target is well inside the m-ball." Therefore, we have: 



In this appendix, a kind of spherical integral, used in the derivation of Equation (30), is given by [38]. 
Let m be a positive integer and / be a real valued measurable function defined on [0, oo), such that 
f(a) > 0 for all a. Then, we have: 



for 0 < r < +oo, where ||x|| is the Euclidean norm on TZ m , and B m is the volume of the unit ball 
in TZ m . 
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